High-fidelity CZ gate for resonator-based superconducting quantum computers 



m 
o 

(N 
£ 

00 



Oh. 



Joydip Ghosh, 1 'B Andrei Galiautdinov, 1, 2 Zhongyuan Zhou, 1 
Alexander N. Korotkov, 2 John M. Martinis, 3 and Michael R. Geller 1 '!]] 

1 Department of Physics and Astronomy, University of Georgia, Athens, Georgia 30602, USA 
2 Department of Electrical Engineering, University of California, Riverside, California 92521, USA 
3 Department of Physics, University of California, Santa Barbara, California 93106, USA 

(Dated: January 10, 2013) 

A possible building block for a scalable quantum computer has recently been demonstrated [M. 
Mariantoni et al., Science 334, 61 (2011)]. This architecture consists of superconducting qubits 
capacitively coupled both to individual memory resonators as well as a common bus. In this work 
we study a natural primitive entangling gate for this and related resonator-based architectures, which 
consists of a controlled-!? 2 (CZ) operation between a qubit and the bus. The CZ gate is implemented 
with the aid of the non-computational qubit |2) state [F. W. Strauch, et al., Phys. Rev. Lett. 91, 
167005 (2003)]. Assuming phase or transmon qubits with 300 MHz anharmonicity, we show that 
by using only low frequency qubit-bias control it is possible to implement the qubit-bus CZ gate 
with 99.9% (99.99%) fidelity in about 17 ns (23 ns) with a realistic two-parameter pulse profile, plus 
two auxiliary z rotations. The fidelity measure we refer to here is a state-averaged intrinsic process 
fidelity, which does not include any effects of noise or decoherence. These results apply to a multi- 
qubit device that includes strongly coupled memory resonators. We investigate the performance 
of the qubit-bus CZ gate as a function of qubit anharmonicity, indentify the dominant intrinsic 
error mechanism and derive an associated fidelity estimator, quantify the pulse shape sensitivity 
and precision requirements, simulate qubit-qubit CZ gates that are mediated by the bus resonator, 
and also attempt a global optimization of system parameters including resonator frequencies and 
couplings. Our results are relevant for a wide range of superconducting hardware designs that 
incorporate resonators and suggest that it should be possible to demonstrate a 99.9% CZ gate with 
existing transmon qubits, which would constitute an important step towards the development of an 
error-corrected superconducting quantum computer. 

PACS numbers: 03.67.Lx, 85.25.Cp 



> 



o 



I. QVN ARCHITECTURE 

Reaching the fidelity threshold for fault-tolerant quan- 
tum computation with superconducting electrical circuits 
[3-[| will probably require improvement in three areas: 
qubit coherence, readout, and qubit-qubit coupling tun- 
ability. Fortunately, the coherence times of supercon- 
ducting transmon qubits 0, [E[ have increased dramat- 
ically, exceeding 10/zs in the three-dimensional version 
d, 0|. Fast, threshold- fidelity nondestructive measure- 
ment has not yet been demonstrated, but is being actively 
pursued [8l-ll2l] . Some method for turning off the interac- 
tion between device elements — beyond simple frequency 
detuning — is also desirable for high-fidelity operations. 
A variety of tunable coupling circuits have been demon- 
strated |13l - fl9{ , but these considerably increase the com- 
plexity of the hardware and it is not clear whether they 
will be practical for large scale implementation. The cou- 
pling can also be controlled by the application of mi- 
crowave pulses [2(il - [32l |. 

An alternative approach has been introduced by 
Mariantoni et al. [33] and theoretically analyzed in 
Ref. [34|. In this quantum von Neumann (QVN) archi- 
tecture, superconducting qubits are capacitively coupled 
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both to individual memory resonators as well as a com- 
mon bus, as illustrated in Fig. [TJ The crossed boxes 
in Fig. Q] represent the phase qubits (35[ employed by 
Mariantoni et al. [33{, however other qubit designs such 
as the transmon may be used here as well. The key fea- 
ture of this architecture is that information (data) is 
stored in memory resonators that are isolated by two 
detuned coupling steps from the bus. Qubits are used 
to transfer information to and from the bus or entangle 
with it, and to implement single-qubit operations, but 
are otherwise kept in their ground states. No more than 
one qubit (attached to the same bus) is to be occupied 
at any time. Such an approach significantly improves 
the effective on/off ratio without introducing the added 
complexity of nonlinear tunable coupling circuitry. The 
spectral crowding problem of the usual qubit-bus archi- 
tecture is greatly reduced because the four-step coupling 
between memory resonators is negligible. And an added 
benefit of the QVN approach is that the longer coherence 
times of the memory elements reduce the overall decoher- 
ence rate of the device. (In Ref. [34[, the architecture we 
consider is referred to as the resonator-zero-qubit archi- 
tecture, but here we will follow the QVN terminology of 
33].) 

The QVN architecture of Mariantoni et al. 33] is not, 
by itself, capable of large-scale, fault-tolerant quantum 
computation, nor is it known how multiple QVN devices 
might be integrated into a scalable design. The problem 
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FIG. 1. Layout of the four-qubit QVN processor. The qi rep- 
resent superconducting qubits capacitively coupled to mem- 
ory resonators m% as well as a resonator bus b. 
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of designing scalable, fault-tolerant architectures for su- 
perconducting qubits is of great interest and importance 
[IS H3, but is still in its infancy. We expect the gate 
design approach discussed here to be applicable to fu- 
ture architecture designs incorporating qubits coupled to 
resonators, and perhaps beyond. 

Along with high-fidelity single-qubit rotations [HI, 
[39j . quantum computing with the QVN processor re- 
quires two additional types of operations. The first 
is state transfer between the different physical compo- 
nents, which has to be performed frequently during a 
computation. The simplest case of state transfer is be- 
tween a qubit and its associated memory (or the reverse) . 
This case is investigated in Ref. [Hj], where two im- 
portant observations are made: First, in contrast with 
the usual SWAP or iSWAP operation, which must be 
able to simultaneously transfer quantum information in 
two directions, only unidirectional state transfers are re- 
quired in the QVN system. This is because adjacent 
qubits and resonators are — by agreement — never simul- 
taneously populated. Second, the phase of a transferred 
|1) state is immaterial, as it can always be adjusted by 
future qubit z rotations [S(| El| . These two simplifica- 
tions allow the resulting state transfer operation, called a 
MOVE gate, to be carried out with extremely high intrin- 
sic fidelity — perfectly for a truncated model — with a sim- 
ple four-parameter pulse profile. By intrinsic fidelity we 
mean the process fidelity (defined below) in the absence 
of noise or decoherence. The need for four control pa- 
rameters immediately follows from the requirement that 
after a MOVE gate, the probability amplitudes must van- 
ish on two device components, the component (q or m) 
the state is leaving, and the bus b. Each zero imposes two 
real parameters, and no other probability amplitudes ac- 
quire weight (in the truncated model). Fixing the phase 
of the MOVE gate, if desired, requires one additional con- 
trol parameter in the form of a local z rotational angle. 

State transfer between a qubit and the bus (or the re- 
verse) can be analyzed in the same way, although in this 
case more pulse-shape parameters are required. In an 
n-qubit QVN processor (consisting of n qubits, n mem- 



FIG. 2. Two-parameter CZ pulse profile ([TJ for the case 
of cj on /27r = 6.8 GHz, u) a/2-K — 7.5 GHz, f ramp = 7ns, 
a = 1.24 ns, and t on = 10 ns. The total gate time excluding 
auxiliary z rotations is t ga to = 17 ns. The example shown 
is representative of a 99.9% fidelity gate for a qubit with 
300 MHz anharmonicity. 



ory resonators, and the bus), zero- amplitude conditions 
must be enforced on the additional n — 1 qubits, lead- 
ing to a total of 2{n + 1) pulse parameters, plus one z 
rotation angle. This makes quasi-exact state transfer to 
and from the bus a considerably more challenging op- 
eration. Simpler three-parameter approximate transfers, 
however, can still be implemented with very high fidelity, 
even when the coupling is strong (see below) . 

Quantum computing with the QVN system also re- 
quires a universal SU(4) entangling gate, the most natu- 
ral being a controlled-cr 2 (CZ) operation between a qubit 
and the bus. The CZ gate investigated here makes es- 
sential use of the non-computational qubit |2) state and 
was first proposed by Strauch and coworkers [42]. The 
Strauch gate has been investigated by man y a uthors and 
has been demonstrated in several systems [43M48I ] . 

The present paper extends previous work by consider- 
ing device parameters and pulse shapes appropriate for 
the QVN system, and by optimizing the CZ fidelity in 
a multi-qubit device. We are especially interested in 
whether the absence of an active tunable coupler re- 
sults in any significant limitation on the obtainable fi- 
delity, given a reasonable amount of qubit coherence, and 
whether very high intrinsic fidelities can be reached with 
a simple and experimentally realistic (filtered rectangu- 
lar) pulse shape. We also study how the gate performance 
rapidly improves with increasing qubit anharmonicity, 
show that the dominant intrinsic error is caused by a 
nonadiabatic excitation of the bus |2) state that occurs 
during the switching the qubit frequency, derive a fidelity 
estimator based on that mechanism, analyze pulse shape 
errors, and simulate qubit-qubit CZ gates mediated by 
the bus. Finally, we address the interesting problem of 
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system optimization, by using gate and idling error esti- 
mates to deduce optimal values of resonator frequencies 
and couplings. 



II. SUMMARY OF RESULTS 

We find that very high intrinsic fidelites — in the range 
of 99.9% to 99.99% and with corresponding total gate 
times in the range of 17 to 23 ns — can indeed be obtained 
with a four-parameter gate. Two control parameters are 



pulse-shape parameters and two are auxiliary local z ro- 
tation angles. We emphasize that only low-frequency 
pulses are required, and that the number of control pa- 
rameters does not depend on the number of qubits in 
the QVN device. The results quoted above assume four 
phase or transmon qubits with 300 MHz anharmonicity; 
other values of anharmonicity are considered below. The 
CZ gate referred to here is between qubit gi and the bus 
(see Fig. [I]), not between a pair of qubits as is usually 
considered. 

The two-parameter low frequency pulse profile we use 
throughout this work is 



€(t) = UJ oK H 
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an example of which is shown in Fig. [5J Here e is the qubit frequency, lo {{ and uj on are the frequencies off and near 
resonance (with the bus), and the pulse switching time is determined by er, the standard deviation of the Gaussians 
inside (fl}. The value of f ramp determines how the pulse is truncated at t = and i ga t as explained in Sec. IIII Fl 
throughout this work we assume that 

Up = 4\/2(j. (2) 

The relation ^ allows the switching time to be alternatively characterized by i ramp , which, as Fig. [^illustrates, is a 
measure of the width of the ramps. The variable i ga tc is the total execution time of the gate excluding z rotations. 
The two control parameters tu on and 

^on — ^-gatc ^ramp (3) 

are determined by the numerical optimization procedure described in Sec. IIII Gl From ([3]) we infer that t on is the time 
interval between the midpoints of the ramps, or the full-width at half-maximum (FWHM) of the pulse. We note that 
the optimal values of t on are somewhat longer than the value 

#suddcn ^ /a \ 

°" ~ V2 9b (4) 

that applies in the sudden, a — ¥ limit. In addition to uj on and t on , two auxiliary local z rotations — on the qubit and 
resonator — are used to implement the CZ gate. As we explain below, adjusting the two control parameters uj on and 
t on zeros the population left in the non-computational qubit |2) state after the gate and (along with the auxiliary z 
rotations) sets the controlled phase equal to —1. The pulse shape JT]) describes a rectangular current or voltage pulse 
sent to the qubit frequency bias through a Gaussian filter of width a, and is believed to be an accurate (although not 
exact) representation of the actual profile seen by the qubits in Ref. [Hj|. 



TABLE I. Optimal state-averaged process fidelity Fa, vo for the Strauch CZ gate between qubit qi and the bus, in the QVN 4 
processor of Fig. [T] No decoherence or noise is included here. Specifications for 99.9% and 99.99% gates are provided for 
three values of qubit anharmonicity n. The parameters i ra m P and a characterize the pulse switching time, and t ga tc is the total 
gate time excluding auxiliary z rotations. -Flu) i s the- minimum fidelity. Data after double vertical lines give the nonadiabatic 
switching error and minimum fidelity estimates; these quantities are defined and discussed in Sec. IIII HI 
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Our main results are given in Table HI Here r\ is the qubit anharmonicity. The 200 MHz results apply to the 
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phase qubits of Ref. [33jj , while the larger anharmonic- 
ities might be relevant for future implementations with 
transmons. The bus couplings g\, are determined by the 
"g optimization" procedure described in Sec. [Vj which 
leads to the simple formula 

— =0.15, (5) 
V 

for the (approximately) optimal bus coupling. Let 
"QVN ra " refer to a quantum von Neumann processor with 
n qubits coupled to n memory resonators and a bus; the 
Hamiltonian for such a device is discussed in Sec. IHI Al 
As indicated in Table U the memory resonators are al- 
ways strongly coupled to allow for fast (less than 5 ns) 
MOVE operations in and out of memory. CZ fidelities 
well above 99.99% are also obtainable (see Sec. IIII J[) . 
Table [J shows that the time required for a qubit-bus CZ 
gate with fixed intrinsic fidelity is inversely proportional 
to the qubit anharmonicity, namely 

(99.9%) 5.2 (99.99%) 6.7 

^ ato sat ° ~V/^' W 

These expressions indicate that CZ gates with very high 
intrinsic fidelity can be implemented in about 20 ns with 
existing superconducting qubits, a conclusion which ap- 
plies not only to QVN„ but also to a wide range of 
similar resonator-based architectures. The intrinsic gate 
(or process) fidelity F avc is the squared overlap of ideal 
and realized final states, averaged over initial states (see 
Sec. IIII E]) . By intrinsic we mean that noise and decoher- 
ence are not included in the gate simulation. The fidelity 
estimate is developed in Sec. IIII HI The results given in 
Table Q] apply specifically to the n = 4 processor, but 
similar results are expected for other (not too large) val- 
ues of n. Two strategies are critical for obtaining this 
high performance: Separating two control parameters 
in the form of auxiliary local z rotations, and defining 
the computational states in terms of interacting system 
eigenfunctions. These strategies were used in Ref. [34| 
and are discussed in more detail below. The gate fideli- 
ties achievable with a transmon-based QVN device are in 



line with that required for fault-tolerant quantum com- 
putation with topological stabilizer codes [49l - [51j . Qubit 
anharmonicity is an important resource that will help us 
achieve that goal. 

The CZ gate of Table U is between qubit q\ and the bus 
in the QVN 4 device, and begins (typically) with a super- 
position of qubit-bus eigenstates, with the other qubits 
and all memory resonators in their ground states. In 
Sec. II VI we also comment on several extensions and vari- 
ations of this basic qubit-bus CZ gate: To begin with, the 
same gate but with qubit 54 (which has a different mem- 
ory frequency) is considered in Sec. lIV Al In Sec. lIVBl wc 
simulate a CZ gate between two qubits in QVN4, starting 
in the idling configuration where the qubits are empty 
and all data is stored in memory. In this case the qubit- 
bus CZ gate of TableQ]is supplemented with MOVE gates 
to effect a CZ between qubits. And in Sec. HVCl we dis- 
cuss the CZ implemented between a pair of directly cou- 
pled anharmonic qubits, instead of a qubit and resonator. 
This is the system originally considered by Strauch et 
al. 0. 

III. CZ GATE DESIGN 

In this section we discuss the qubit-bus CZ gate design 
problem. 

A. QVN model 

The QVN n processor consists of n superconducting 
qubits [H4H capacitively coupled [42],[52|,[53| to n memory 
resonators and to a common bus resonator. Here we as- 
sume parameters appropriate either for phase qubits (35j | 
or transmon qubits 4, 5] with tunable transition frequen- 
cies. We write the qubit angular frequencies as ej, with 
i = 1, • • • ,n. These are the only controllable parameters 
in the QVN Hamiltonian (in contrast with Refs. [I3-[32j 
we do not make use of microwave pulses). The memory 
frequencies are written as w m j, and the bus frequency is 
cob- The (bare) frequencies of all resonators are assumed 
here to be fixed. 



Because we are interested in very high fidelities, a realistic model is required. However, we have shown (in un- 
published work) that the CZ performance is extremely robust with respect to the model details, so we only report 
results for a simplified Hamiltonian; the approximations used are discussed below. For the qubit-bus CZ simulations, 
the Hilbert space is truncated to allow for up to three excitations. The CZ gate naively involves no more than two 
excitations, so to properly account for leakage we include up to three. Therefore, four-level qubits and resonators 
(which include the |3) states) are required in the model. The QVN Hamiltonian we use is 
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excluding single-qubit terms for microwave pulses that 
are not used in this work. Here r\ and 77' are qubit an- 
harmonic detuning frequencies, g m and gb are the qubit- 
memory and qubit-bus interation strengths, and 
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The matrices in (JJJ act nontrivially in the spaces indi- 
cated by their subscripts, and as the identity otherwise. 
The matrix Y results from a harmonic oscillator approx- 
imation for the qubit eigenf unctions. Factors of h are 
suppressed throughout this paper. 

The main approximations leading to ([7]) are the neglect 
of the e-dependence of the interaction strengths g m and 
gb, and the neglect of a small direct coupling between the 
memories and bus [54j . We have verified that including 
these does not change the main conclusions of this work. 
The e-dependence of the anharmonicities 77 and 77', and 
small anharmonic corrections to the interaction terms in 
(0, are also neglected. 

The parameter values we use in our simulations are 
provided in Table [TT1 We assume 77' = which is appro- 
priate for qubic anharmonicity. As discussed in Sec. |Vl 
the value of the bus coupling gb is chosen to give the 
shortest CZ gate time (for a range of fidelities). The 
choice of resonator frequencies is also discussed in Sec.fVl 
We simulate n = 4 qubits. The fidelities quoted in 
this paper are numerically exact for the model ([7]); the 
rotating-wave approximation is not used. 



TABLE II. Device parameters used in this work. 



quantity 


value 


empty qubit parking frequency u par k/27r 


10.0 GHz 


memory resonator mi frequency LJmi/2% 
memory resonator 7712 frequency LUm2/2n 
memory resonator 771,3 frequency ui-mz^n 
memory resonator 777.4 frequency cu„ l 4,/2n 


8.3 GHz 
8.2 GHz 
8.1GHz 
8.0 GHz 


initial detuned qubit frequency ui g/2n 


7.5 GHz 


bus resonator frequency uj^/2-k 


6.5 GHz 


qubit-memory coupling strength g m /2ir 


100 MHz 


qubit-bus coupling strength gb/2ir 


30 - 60 MHz 


qubit anharmonicity 77/27r 


200 - 400 MHz 



Although the CZ and MOVE gates considered here do 
not involve microwave pulses, the single-qubit gates are 
assumed to be implemented with microwaves in the usual 
manner at the qubit frequency Lu g. This frequency is 
also used to define an experimental "rotating" reference 
frame or local clock for each qubit: All qubit frequencies 
are defined relative to ui s [33[. This is discussed below 
in Sec. HITCl 



B. Strauch CZ gate 

In this section we give a detailed description of the CZ 
gate introduced by Strauch et al. [42J. In particular, we 
explain the specific roles played by the two pulse-shape 
parameters t on and uj on , and by the two auxiliary z rota- 
tion angles 71 and 72. To accomplish this we introduce 
several approximations that allow for an analytic treat- 
ment of the CZ gate dynamics. 

First, we consider a truncated model consisting of a 
single superconducting qubit with frequency e and an- 
harmonic detuning 77, capacitively coupled to a bus res- 
onator with frequency u>b, 




H 



In this case Y reduces to 
'0 



gbY q (g>Y h . (9) 
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This Hamiltonian is written in the basis of bare eigen- 
states, which are the system eigenfunctions when the 
qubit and resonator are uncoupled. We write these bare 
states as \qb), with q,b £ {0,1,2}. The energies of the 
interacting eigenstates, which we write with an overline 
as I qb), are plotted in Fig. [3] as a function of e/2ir for 
the case of uj h /2-K = 6.5 GHz, 77/2^ = 300 MHz, and 
(?b/27T = 45 MHz. The interacting eigenstates are labeled 
such that \qb) is perturbatively connected to \qb) when 
e > 0Jb- 

Second, we assume a short switching time and ignore 
the dynamical phases acquired during the ramps. As we 
will see below, this approximation is valid when gb -C 77, 
so that the switching can be made sudden with respect 
to the coupling gb, but still adiabatic with respect to the 
anharmonicity 77. 

The CZ gate of Strauch et al. |42j |. adapted to the qubit- 
resonator system, works by using the anticrossing of the 
1 11) channel with the auxiliary state |20). In terms of 
the pulse parameters defined in ([1]), the qubit- resonator 
state is prepared at a qubit frequency e = L) g, and the 
frequency is then switched to e = u on for a FWHM time 
duration 7j on . In the simplified model considered in this 
section, 



and 



Won = W b + 77, 



tr 



V2~9h 



(11) 



(12) 



Equation (jlll) gives the qubit frequency for which the 
bare state |11) is degenerate with |20), and is at a fre- 
quency 77 above the usual resonance condition. Equa- 
tion (fl"2"|) is the sudden-limit value defined in (Q| and de- 
rived below. The qubit frequency is then returned to 



6 



N 
X 

O 

•« 

E» 

<D 
c 



14.8 
14.4 
14.0 
13.6 
13.2 
12.8 
12.4 
I2.0 
7.6 
7.2 
6.8 
6.4 
6.0 





|20) 






^ 




02} 








.s ■ 

\0l) 








foo) 



6.50 6.75 7.00 7.25 
qubit frequency (GHz) 



7.50 



FIG. 3. (Color online) Energies of eigenstates \qb) of a single 
qubit q coupled to a resonator bus b. Here u;t>/27r = 6.5 GHz, 
■q/2-K = 300 MHz, and g h /2n = 45 MHz. The time depen- 
dence of the qubit frequency during a CZ gate (solid black 
curve) is indicated at the top of the figure. 



the detuned value w ff- The complete pulse profile is 
also shown in Fig. [3] (solid black curve) for the case of 
ijj ori /2Ti — 6.8 GHz and w ff/27r = 7.5 GHz. 

Let's follow the evolution resulting from an initial (nor- 
malized) qubit-resonator state 



|00> + a i|01> + aio|l0>+aii|ll>. 



(13) 



Because the 1 00) channel is very well separated from the 
others, the |00) component will only acquire a dynamical 
phase factor 



(14) 



where Eqq is the energy of the |00) eigenstate. Without 
loss of generality we can shift the entire spectrum so that 
E 00 = [as in ©] and the phase factor (IT4"|) becomes 
unity. This freedom results from the fact that any unitary 
gate operation only needs to be defined up to an overall 
multiplicative phase factor. With this phase convention 
the CZ gate acts as the identity on this component, so 
we have the map 



The 1 01) component will mostly return to |01), also 
with an acquired phase, but a small component will be 
left in 1 10) due to the nonadiabatic excitation of that 
channel, which is only separated in energy from |01) by 
about rj when e = u>h + V- The 1 10) component similarly 
suffers from a small noniadabatic coupling to |01). As 
we will explain below, these nonadiabatic errors are ex- 
ponentially suppressed when the functional form of e(t) 
is properly designed. Then we have 



|01) 



Vl-IEi |01) +e- ia 'y/E~ 1 1 10) (16) 



and 



1 10) 



-id 



V 7 ! - Ei 1 10) + e^'-v/Ei |01), (17) 



where Ei is a small nonadiabatic population error (below 
we refer to Ei as a switching error). In the Ei — > limit, 
a and j3 are dynamical phases given by 
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E 01 dt : 



E w dt : 



9l 



(18) 
(19) 



where the second approximate quantities neglect phase 
accumulation during the ramps and use perturbative ex- 
pressions for the energies Eqi and E\q when e = ujb + V- 
The expressions (fT6|) and (flT|) neglect an extremely small 
leakage out of the {|01), |10)} subspace. Neglecting this 
leakage, the evolution in the {|01), |10)} subspace is uni- 
tary, leading to the phase condition 



0. 



(20) 



Using (p?0|) to eliminate j3' leads to 

Vl-Ei jbT) + e ^+^ ^17 JT0>, (21) 



01) 



1 00) ->■ 1 00). 



(15) 



1 10) ->■ e _i Vl - Ei 1 10) - e-^-^V^ET |01>, (22) 

where <j> = a' — (3. The evolution of the eigentates 1 01) 
and 1 10) is therefore characterized by the cross-excitation 
probability Ei and three phase angles a, /?, and <p. 

Now we consider the |11) component. The |11) channel 
couples strongly with the 1 20) channel, as well as weakly 
with 1 02). The simplest way to understand the dynamics 
of the 1 11) component is to use two different represen- 
tations to describe thes e tw o effects. We will describe 
strong interaction with 1 20) in the bare basis and the 
weak, nonadiabatic coupling with 1 02) in the eigenstate 
basis. Suppose we be gin w ith the qubit strongly detuned 
from the bus, so that 1 1 1) ~ 1 1 1 ) (the detuned interacting 
eigenstate is well approximated by the bare |11) state). 
Then we quickly switch e from tu g to Wb+?7- By "quickly" 
we mean that we strongly mix with the 1 20) channel. The 
interaction with 1 02) is always weak, even in the sudden 
limit. This asymmetric excitation is possible because 1 20) 
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is protected (separated in energy from 1 11)) by an energy 
gap 2v / 2gb, whereas |02) is protected by a much larger 
gap of 77 — V2gb (this expression accounts for level repul- 
sion from 1 20), and we have assumed that g b <C r/). We 
can informally say that the desired switching is nonadia- 
batic with respect to the energy scale g b , but is adiabatic 
with respect to rj [42J. 

Focusing hrst on the strong coupling to |20), the sud- 
denly switched |11) state is no longer an eigenstate when 
e = uj oni as the relevant eigenfunctions at this setting are 

|rT)J n) - |20) and P)= |11) + ' 20) . (23) 



V2 

The nonstationary state 



V2 



UD = 



111) 



|20) 



V2 



(24) 



therefore rotates in the {|11), 1 20) } subspace, and after a 
time duration t becomes 



-iEnt 



-iEnt 



111) 



iAEt 



|20) 



V2 
1 + e- lAEt 



1 



-iAEt 



1 20) 



where 



AE = E 2Q - E u = 2V2g b . 



(25) 



(26) 



(27) 



Holding e fixed at uj b + rj for a F WHM time (fT2~j) , corre- 
sponding to a 2ir rotation, (j2"6")l becomes 



-i-Ellton 



11) 



When e = Ldb + rj, the energy of eigenstate |11) is 

En = 2uj h + rj- V2g h . 
After detuning quickly we therefore obtain 

_2wb + rf 



|H) 



exp 



V2g b 

or, using expressions (fTS)) and (H"9l) . 



|H). 



|11) -> -e-^+V |n). 



(28) 



(29) 



(30) 



(31) 



The two phase angles a and /3 can be cancelled by 
the application of independent auxiliary single-qubit z 
rotations 



R z (j) = exp[- l ( 7 /2K] 



(32) 



to the qubit and bus. Qubit z rotations are implemented 
by frequency excursions, whereras resonator z rotations 



are implemented in software (they are compiled into fu- 
ture qubit rotations). Following the pulse sequence that 
leads to ([T5]). (j2T]t . ([22]). and ([3T]>. with the operation 



where 



fl,(7i)8>-R*(*7!2), 



-/3 and 72 = —a, 



leads to the map 



00) ->■ |00), 



01) -> Vl-Ei 1 01) + e-^VEi 1 10), 
10) ->■ Vl-Ei jlO) - e^-s/Wx ]oi), 



|11)^-|11), 



(33) 
(34) 

(35) 
(36) 
(37) 
(38) 



apart from a global phase factor. The use of auxiliary z 
rotations is discussed further Sec. IIIIDI 

The minus sign in (|3"5)) is the key to the Strauch CZ 
gate. However, as mentioned above, the analysis leading 
to (f3"8"| neglected a weak nonadiabatic excitation of the 
1 02) channel caused by the switching of e. Including this 
effect in (l38l) leads to the modification 



|11> Vl-Ea |H) 



phase factor x VE2 |02), 



(39) 



where E2 is another switching error. Both Ei and E2 van- 
ish exponentially with a (or £ ra mp), and for the regimes 
studied in this work E2 is the dominant source of intrin- 
sic gate fidelity loss. We note that the analysis leading 
to (|3"9")) assumed implementation of the ideal values [ ([TTj) 
and (|12[) ] of w on and t on . Errors in these two control 
parameters, which we refer to as pulse shape errors and 
study in Sec. IIII II lead instead to 



|ll) ->■ - e'Vl -E 2 -K g |ll) 



phase factor x \/E2 02) 
phase factor x yEe |20), 



(40) 



where the controlled-phase error angle (5 and rotation er- 
ror Eg depend on the errors in w on and t on , respectively. 

Finally, it is also interesting to consider the fully adi- 
abatic limit of the Strauch CZ gate. By this we mean 
that the switching is adiabatic with respect to both g b 
and rj. For the gate time to be competitive with the 
nonadiabatic gate of Table HI a larger coupling <?b is re- 
quired, which might lead to significant higher-order and 
cross-coupling errors in a multi-qubit device, but in the 
fully adiabatic limit only one pulse control parameter — 
either ui on or t on — needs to be optimized (two z rotations 
are still required). This is because adiabaticity now as- 
sures that the |11) population is preserved (apart from 
exponentially small switching errors), taking over the role 
previously played by i on , and a single pulse shape param- 
eter is sufficient to specify the controlled phase. A highly 
adiabatic CZ gate was demonstrated in Ref. |43| . 
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There are a few important differences between the 
Strauch CZ gate applied to a pair of directly coupled 
qubits (as in Ref. [421 ]) and to the qubit-bus system con- 
sidered here. These differences result from the harmonic 
spectrum of the resonator in the latter case and are dis- 
cussed below in Sec. IIV CI 

C. Eigenstate basis 

The Hamiltonian (|7|) is written in the usual bare basis 
of uncoupled system eigenstates, but information pro- 
cessing itself is best performed in the basis of interacting 
eigenfunctions of H idlc , where H[ d \ e is given b y Q with 
the qubits in a dispersive idling configuration [34]. This 
choice of computational basis assures that idling qubits 
suffer no population change in the decoherence-free limit, 
and evolve in phase in a way that can be almost exactly 
compensated for by an appropriate choice of only n rotat- 
ing frames or local clocks, one for each qubit [34]]. Here 
we briefly review this important concept. 

In principal, any complete orthonormal basis of the 
physical Hilbert space that can be appropriately pre- 
pared, unitarily transformed, and measured — essentially, 
any basis where one can implement the DiVincenzo crite- 
ria [55l ] — is a valid basis on which to run a quantum com- 
putation. Defining the computational states to be inter- 
acting system eigenfunctions gives them the simplifying 
property that the time evolution can be decomposed into 
a sequence of gates, between which (almost) no evolution 
occurs. In other words, idling between gates generates 
the identity operation. This property, which is implic- 
itly assumed in the standard circuit model of quantum 
computation, could be realized in an architecture where 
the Hamiltonian H can be completely switched off be- 
tween gates. However, it is not possible to set H = in 
the QVN architecture; nor can H itself be made negligi- 
bly small between gates. Therefore, nonstationary states 
such as uncoupled-qubit eigenstates accumulate errors 
(including population oscillations) between gates unless 
a correction protocol such as dynamical decoupling [56] 
is used. By defining computation al st ates in terms of 
interacting system eigenfunctions a t some prede- 

fined dispersive idling configuration (qubit frequencies), 
the only evolution occuring during an idle from time t\ 
to t 2 is a pure phase evolution, 

m^^i^^^^^Wh)), (4i) 

where E is the exact energy eigenvalue (and we neglect 
dccohcrcncc). Furthermore, it is possible to compensate 
for — or effectively remove — the pure phase evolution in 
(HH) by applying phase shifts (after the idle period) to 
each eigenfunction to cancel the e -lB ( t2-tl ) phase factors; 
doing so would result in the ideal between-gate evolution 

W^^WM^Wm))- (42) 

The idling dynamics ([42]) is evidently equivalent to set- 
ting H = between gates. We will discuss below how 



the compensating phase shifts are actually implemented 
in practice. 

This use of interacting system eigenfunctions and com- 
pensating phase shifts as described above provides a com- 
putational basis that evolves ideally between gates, but 
such an approach is not scalable; for example, there are 
2 2n+1 such computational states in QVN n . In Ref. [H| 
an approximate but scalable implementation of this ap- 
proach was introduced. The idea is that the exact energy 
E of a computational state in QVN„ is, to an extremely 
good approximation, the sum of uncoupled qubit and res- 
onator frequencies, i.e., essentially noninteracting. This 
is not simply a consequence of the dispersive regime en- 
ergies (eigenvalues of i/idic), which have non-negligible 
interaction corrections, but because only a special subset 
of the eigenfunctions are used for information processing: 
In the QVN„ system we only make use of -ffidio eigen- 
functions in which there are no more than n excitations 
present, and such that two directly coupled elements — 
qubits or resonators — are not simultaneously occupied 
(except during the CZ gate). For example, when the 
data is stored in memory, the residual memory-memory 
coupling is fourth order in the qubit-resonator coupling g 
(for simplicity we assume here that <?b = 3m)- This leads 
to an eighth-order conditional frequency shift (order g 16 
idling error) (34{. Next, suppose an excitation is trans- 
ferred from memory to a qubit via a MOVE gate. Now 
the dominant frequency shift is sixth order. And when 
an excitation is in the bus the largest shift is fourth order 
[34] ]. The largest idling error (associated with the phase 
compensation) is therefore eighth order in g and can be 
made negligible with proper system design. 

The compensating phase shifts could be implemented 
through additional local z rotations, one for each qubit 
and resonator. However, these phase shifts evolve in 
time with very high (> I GHz) frequency, and it is there- 
fore experimentally more practical to introduce a local 
clock/rotating frame for each qubit and resonator. This 
is achieved by introducing a fixed-frequency microwave 
line for each qubit and resonator, and measuring each 
qubit and resonator phase relative to the phase of its 
reference. By choosing the frequency of the qubit (res- 
onator) reference microwave equal to the idle frequency 
(resonator frequency), the component frequencies [and 
therefore the quantity E in (|4"1"T) ] are effectively zeroed, 
and no more than 2n + I different reference frequencies 
or local clocks are required. This procedure corresponds 
to implementing the experiment in a multi-qubit rotat- 
ing frame. And, in a further simplification, the local 
clocks/rotating frames for the resonators are replaced by 
additional qubit z rotations that are handled in software 
(i.e., combined with future rotations). Therefore, in prac- 
tice only n local clocks/rotating frames are needed, one 
for each qubit. 

Because the CZ gate simulations reported in Table U 
are already supplemented with local z rotations, these 
local clocks/rotating frames do not need to be included 
in those simulations; we simulate the lab frame. How- 
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ever, they are included in the pulse-shape error simula- 
tions reported after (|9~4")l and in the CZ23 qubit-qubit gate 
simulation reported in Sec. IIVBI 

Having motivated the use of interacting system eigen- 
functions for computational basis states, it is still neces- 
sary to establish that such states can actually be prepared 
and measured. Because we can assume the processor to 
initially start in its interacting ground state — a compu- 
tational basis state — preparation of the other computa- 
tional states can be viewed as a series of n pulses and 
MOVE gates. We expect that such operations on the 
interacting eigenfunctions can be performed at least as 
accurately as when applied to bare states. Eigenfunc- 
tion readout is a more subtle (and mo del- dependent) 
question, but the analysis of Ref. 34] suggests that 
interacting-eigenfunction readout is actually better than 
bare-state readout (in the model considered there). 

We also note that the idling configuration and associ- 
ated eigenstate basis generally changes between consec- 
utive gates (an example is given below in Sec. IIVBI) . In 
Table HI the idling configuration has qubit q\ at w ff and 
the others at w par k. Therefore, our entangling gate de- 
sign is constrained by the requirement that we start and 
end in eigenstates of this particular -ffidie- 

The discussion above motivating the use of interacting 
eigenstates is based on their nearly ideal idling dynam- 
ics. It is still interesting, then, to consider whether the 
CZ gate can be generated equally well in either (bare 
or interacting eigenfunction) basis. We find that for the 
parameter regimes considered here, it is not possble to 
achieve better than about 99% fidelity in the bare ba- 
sis with the same two-parameter pulse profile (it should 
be possible using more complex pulse shapes). The re- 
maining error is consistent with the size of the perturba- 
tive corrections to the bare states in the idling configu- 
ration. This exercise emphasizes the importance of per- 
forming quantum logic with the system eigenfunctions, 
which have the built-in protection of adiabiticity against 
unwanted transitions. 

One might object to the use of interacting eigenfunc- 
tions as a design tool, the exact calculation of which is not 
scalable. However, approximate dispersive-regime eigen- 
functions are efficiently computable. A particularly sim- 
ple way to do this is to calculate the generator S of the 



in 9b and g m . At the 99.99% fidelity level, it is sufficient 
to calculate S to first order. Writing H i( u c = H + 5H 
leads to the condition i[S, Hq] + 8H = 0, which is im- 
mediately solvable in the bare basis \q1q2 • • ■ TO1TO2 •••&). 
Here qi,mi,b € {0, 1, 2, ... }. Other efficient eigenfunc- 
tion approximation schemes are also possible. 

In this work we denote the exact or approximate H-^ic 
eigenfunction perturbatively connected to the bare state 
\q1q2 ■ ■ ■ q n mim 2 ■ ■ ■ m n b) by 



\q1q2 ■ ■ ■ qnm 1 m 2 ■ ■ ■ m n b) , 



(43) 



following the overline notation introduced above. Note 
that (|43l) is not a tensor product of single-qubit / resonator 
eigenstates as is usually the case. 



D. Auxiliary z rotations and CZ equivalence class 

The standard CZ gate in the bare two-qubit basis 
{|00),|01),|10),|ll)}is 



CZ: 



a 
0100 
0010 
\0 -1, 



(44) 



However in the QVN processor, local z rotations can be 
performed quickly and accurately, typically by brief qubit 
frequency excursions. Thus, we will consider the limit 
where SU(2) operations of the form exp[— i(6/2)a z ] can 
be done on the qubits and bus with negligible error and 
in a negligible amount of time (fidelity loss resulting from 
errors in these rotations are discussed in Sec. IHIlj) . We 
therefore want to define our entangling gate modulo these 
z rotations. We will do this by constructing a local-z 
equivalence class for an arbitrary element (gate) in SU(4), 
and then specialize to the CZ gate. 

We define two elements U and U' of SU(4) to be equiv- 
alent, and write U' = U, if 



U' 



A post 



UUr; 



(45) 



diagonalizing transformation V 



-iS 



by a power series where 



U ( 71 , 72 ) = R Z ( 71 ) <g> R Z { 12 ) = e ^+^)/ 2 



(1 



























e -*n 













e -i(7l+72)y 



(46) 



for some rotation angles 7^. The local-z equivalence class {U} corresponding to U is the set of elements u p0 st U u pYC 
for all itprc, Upost- For a given gate U, {U} typically occupies a four-dimensional manifold, depending on four rotation 
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angles. But because (|44|) is diagonal, {CZ} instead forms a two-dimensional sheet, 



{CZ} = phase factor x 



/ 1 

e- l ~< 2 



\0 







-4(71+72) 



(47) 



The CZ gate (|44[) can be obtained by reac hing any point in the {CZ} plane and then performing auxiliary z rotations. 
And it is straightforward to confirm that [43 



(-1 0\ 

10 

10 

V ooi/ 




TO 0\ 
10 
0-10 

v0 1/ 




(48) 



We note that bus rotations, which cannot be directly 
implemented with microwave pulses or frequency excur- 
sions, are compiled into future qubit rotations. 

The discussion above assumed a pair of qubits or a 
qubit and resonator, but it applies to a QVN processor 
in the interacting eigenfunction basis (I43[) after a minor 
modification. In the bare basis, the CZ gate is typically 
defined through its action (l44l) on a pair of qubits (or a 
qubit and resonator). Then, action on a bare computa- 
tional basis state such as \q1q2 • • • (Irjnvmi • ■ ■ m n b) fol- 
lows from the tensor-product form of that bare state. In 
the eigenstate basis the CZ gate must be defined through 
its action on 



\qiq2 ■ ■ ■q n mim 2 ■ ■■m n b), 



(49) 



such as to reproduce the ideal action on the bare states 
to which they are perturbativey connected. For example, 
the CZ gate on qubit qi and the bus acts ideally as 



CZ \Oq2q3q4m1m2m3m4O) = \Oq2q3q4m1m2m3m4O) 



CZ \Qq2q3q4m1m2m3m4l) = \Qq2q3q4Tn1m2rn3rn4X) 



CZ \ lq2q3q4m1m2m3m4Q) = \lq2q3q4m1m2m3m4O) 



CZ \ \q2q3q4m\m2m3m4\) = - |lg2<?394"iim 2 m 3 rn4l), 
where qi,nii € {0, 1,2,...}. 



^target = CZ [see flUJl]. Note that the projected U is 
not necessarily unitary here, and that the first term in 
([50)1 characterizes the possible leakage from the compu- 
tational basis (non-unitarity) whereas the second term is 
proportional to the square of the Hilbert-Schmidt inner 
product of U with [/target ■ Although U is not assumed to 
be unitary, the expression (|50[) assumes a pure state and 
is (obviously) not valid in the presence of decoherence. 
[The formula assumes that the Kraus representation 
for the completely positive process is not necessarily trace 
preserving, but it has only one term.] The form (1501) also 
assumes an average over a 4-dimensional Hilbert space; 
in the iV-dimcnsional generalization the denominator is 
N + N 2 , which is necessary (note numerator) to assure 

that F avc ([/target, C/target) IS Unity. 

It is also useful to calculate the minimum or worst-case 
fidelity. The minimum fidelity of interest here is the state 
fidelity minimized over initial computational states, for 
a gate that has already been optimized (by maximizing 
-Pave). In Sec. IIII Bl we argued that the dominant intrinsic 
error mechanism (for an optimal pulse) in the truncated 
qubit-resonator model ^ is the nonadiabatic excitation 
of the 1 02) channel, in which there are two photons left 
in the bus resonator. Therefore, in the model (|9]), the 
min imum state fidelity occurs for the initial eigenstate 
|11). In the QVN4 processor, this worst-case state is 
written [in the notation of (l43l) ] as 



E. Fidelity definitions 

The gate or process fidelity measure we use in this work 
is based on a state fidelity defined by the inner product of 
the ideal and realized final (pure) states, scmared. This 
leads to a state-averaged fidelity given by [57], [H[ 



-^ave (^> ^target) 



Tr([/t[/) + |Tr([4 rget [/)| 
20 



(50) 



where U is the realized time-evolution operator in the in- 
teracting eigenfunction basis after auxiliary z rotations, 
projected into the relevant computational subspace, and 



1100000001), 



(51) 



where we have assumed a CZ gate between qubit q\ and 
the bus. Numerical simulation of this gate in the QVN4 
processor confirms that the minimum CZ fidelity indeed 
occurs for the initial state (IBTI) . and is due to leakage 
from the computational subspace. We therefore define 
the minimum fidelity to be the state fidelity for initial 
condition (|5"TT). 



F lll} = I (100000001 |t/| 100000001) I ' 



(52) 



Note that this expression is not sensitive to the value of 
the controlled phase, and only accounts for leakage from 
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the computational subspace. Simulated values of -F|n) 
for the qubit-bus CZ gate are given in Table HI along 
with estimates of this same quantity that are discussed 
below in Sec. UTTHl 



F. Pulse shape 

In the QVN Hamiltonian (J7|), the qubit frequencies €i 
are the only available experimental controls. [There are 
also single-qubit terms for microwave pulses that are not 
shown in ((7|) and not used in this work.] During a CZ 
gate between a given qubit and the resonator bus, the 
frequency of that qubit is varied according to ([T]) , where 



Erf(t) 



~ x dx. 



(53) 



with the other qubits remaining at the parking frequency 
(given in Table [n} . 

Two quantities related to the pulse switching — a and 
^ramp — appear in the pulse profile function ((T|). We do 
this to emphasize that, in principal, two independent 
quantities could be used to quantify the shape of the 
switching profile. The first is the time duration (or 
width) of the switching, characterized by the standard 
deviation a. The second is the trunction time of the 
pulse, measured from the center (half maximum) of the 
frequency switch, which is equal to t lamp /2. However in 
this work we always use the (somewhat arbitrary) rela- 
tion ©, which amounts to cutting off the pulse at 2\/2 
standard deviations from the switching midpoints. 

If the pulse shape function (TTJ) is used only when < 
t < igatc, and e is set to uj a g otherwise, there will be 
small pulse discontinuities at t = Q and t gatc , the size 
of which is determined by the condition Assuming 
^ramp < ^gato, we have 



e(0) = e(tg a to) « w off + 



1 - Erf (2) 



(54) 



which differs from the asymptotic value u> g by an 
amount 



Woff 



1 - Erf (2) 



(55) 



where 1 — Erf(2) sa 0.5%. However, in an experiment 
these discontinuities are usually smoothed over by addi- 
tional pulse shaping. Moreover, our simulations begin 
t = and end at t — t gatc , so the truncation only slightly 
affects the pulse shape: The initial and final detuned 
qubit frequency is actually a few MHz smaller than w ff . 

Having fixed the relation @, there are then two mea- 
sures of the pulse switching time, a and i ra mp, with i ramp 
providing a convenient measure of the time duration of 
the ramps. This property can be seen in the pulse shape 
example of Fig. 



G. Gate optimization 

Wc numerically optimize the two pulse-shape control 
parameters w on and t on , as well as the two auxiliary z 
rotation angles, 71 and 72, to maximize the average gate 
fidelity (|50l) . All other pulse parameters (uj g, a, t ramp ) 
are fixed. The values of w g and £ ra mp are determined 
by the system optimization analysis of Sec. El and then 
is obtained through relation @. The roles played by 
the control parameters u; on and t on are discussed above 
in Sec. UTTBl 

The fidelity optimization procedure is carried out in 
two stages: In the first stage we take w on to be equal to 
its approximate value [see (fTTT) ] 



(56) 



and optimize i on to get close to the two-dimensional 
equivalence class {CZ} defined in (|4T|) . We do this by 
minimizing a sum of two positive errors, one measur- 
ing the deviation of the absolute values of the matrix 
elements of the evolution operator U (projected into the 
qi~b subspace) from that of the four-dimensional identity 
matrix, the other measuring the deviation from the ideal 
relationship between the phases of the diagonal elements 
indicated in (|4"T|) . This first stage yields an approximate 
value of t on , as well as approximate rotation angles 



71 «arg (10|C/|10) 

72 « arg W\UW) 



arg (00|[/]00), 
arg (OOjC/fOO). 



(57) 
(58) 



In the second stage of optimization, we use the ap- 
proximate values of uj on , t on , 71, and 72, obtained from 
the first stage, as seeds for a full four-dimensional (uj m 
tom 7ij 72) nonlinear maximization of the average fidelity 
(1511 between 



■"(71,72) x U 



(59) 



and the standard CZ gate (|44[) . Here U is the projected 
evolution operator and u € SU(2) <£> SU(2) is defined in 
(H2). 

There are two obvious approaches to simulating the 
gate evolution in the interacting eigenfunction basis (|43l) . 
The first is to use the appropriate i?idic eigenstates to 
transform the target gate (|4"4"|) to the bare basis, and 
perform the actual simulation in the bare basis. How- 
ever this approach requires that the subspace projection 
operators — which are naturally defined in the eigenstate 
basis — also be transformed. The second approach, which 
we follow here, is to transform the QVN Hamiltonian ([7]) 
to the eigenstate basis and perform the simulation in that 
basis. 



H. Switching error and fidelity estimate 

In this section we calculate the transition probability 
caused by a change of the qubit frequency during a CZ 
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pulse, MOVE gate, or any other operation in the QVN 
system or related superconducting architectures. The 
problem will be treated quite generally and then applied 
to the Strauch CZ gate of Table HI resulting in a simple fi- 
delity estimator for that gate. The quantities A, p sw , and 
appearing in Table Q] are discussed in this section. 

Imagine that we have prepared an initial interacting 
system eigenfunction \a) prior to performing a CZ oper- 
ation or other gate that involves changing the frequency 
of one or more qubits. We assume that the ideal (tar- 
get) behavior during the frequency switch or ramp itself 
is the identity map (times a phase factor), and that the 
| a) channel does not cross any others in the system. The 
population loss during the ramp will therefore be expo- 
nentially suppressed if the switching time is long enough. 

In a multi-qubit system there is typically a large num- 
ber of nonresonant channels coupled to \a) that can be ex- 
cited by the frequency switch. However, when the ramp 
fidelity is high and the probabilities of the undesired tran- 
sitions 



\a)^\b),\b>),\V>),- 



(60) 



are small, they can be individually estimated pertur- 
batively (neglecting interference), thereby reducing the 
problem to a sum of independent two-channel problems 



Jo) 
Jo) 
la) 



\b») 



(61) 



each characterized by a time-dependent detuning A and 
a coupling G. Without loss of generality we can shift the 
energy of a given two-channel problem so that the bare 
final state has zero energy. Each nonadiabatic transi- 
tion can therefore be described by a general two-channel 
model of the form 



H = 



A(t) G 
G 



(62) 



in the bare basis spanned by {|o),|6)}. The undesired 
final state |&) has a fixed energy and the energy of 
| a) varies in time with detuning A. The coupling G is 
assumed to be a real, positive constant. 
The instantaneous eigenstates of (|62[) are 



where 



I a) = cos f 1°) +sin||6), 
1 6) = cos -| 6) - sin-|a), 



,'2G\ 

\ = arctan [ — I . 



The instantaneous energies are 




(63) 
(64) 

(65) 
(66) 



and 




(67) 



The \b) channel is initially unoccupied at time t = 0, and 
we are interested in the probability p sw that the system, 
prepared in |a), is found in \b) after changing the de- 
tuning A from one value to another. We refer to this 
probability as the nonadiabatic switching error, which 
we calculate by expanding the wave function in the basis 
of instantaneous eigenstates (|63[) and (|64l) , as 



m—a.b 



,e-*/o*«".*-| m ). 



This leads to 
dt 

where 



^dA W dt 



77T d TT G 
(b — — o) = --r-T x 
x >dA> 1 A 2 



1 



(68) 



(69) 



(70) 



The nonadiabatic matrix element (|70p has been written 
so that the second term approaches unity in the G< A 
perturbative limit. 

At time t — 0, tp a = 1. An approximate expression for 



IV'ft^final)^ 



(71) 



can be obtained from (|69p by assuming that \ipf,\ <C 1 
throughout the evolution, so that for all t. Then 



Psw — 

where 



GA« 



O 2 



-dt = 



Xe- l ti ndT dt 



fl = E a -E b = VA 2 + 4G 2 



(72) 



(73) 



is the instantaneous splitting. We can simplify (|72p fur- 
ther by assuming G <C A, which will be the case for the 
applications considered below. In this perturbative limit 
we therefore obtain 



Ps 



GA 



e -ifjAdT dt 



(74) 



We emphasize that the form (174)) assumes that A does 
not pass through zero, which would cause Landau-Zener 
tunneling and invalidate the perturbative analysis. 

In this work we are specifically interested in p sw for a 
single switch of the detuning from A on to A g (or the 
reverse) according to the smooth, error-function based 
profile 



A(t) 



A off + A 



A off - A 



■Erf 



f - it 



2 'ramp 



V2cr 



(75) 
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FIG. 4. Detuning pulse profile ([75jl for a single frequency 
switch. 

shown in Fig. [4] The standard deviation a character- 
izes the switching time of the pulse and i ra mp [related 
to a through j2])] specifies its truncation, as discussed in 
Sec. IIIIFl We use this switching profile for both CZ and 



MOVE gates. The switching error for the single switch 
profile defined in (|75|) and shown in Fig. |4] can be ex- 
pressed as 

^=(£)V' (76) 

where 

A(A on ,A off ,a) ee Ao-I^e-*/.'^^! (77) 

The dimensionless quantity \A\ 2 is plotted in Fig. [5] for 
five instances of A on and A D ff relevant to this work. We 
note that p sw evidently decreases as an exponential func- 
tion of a, as expected for a nonadiabatic process. How- 
ever, the dependence of \A\ 2 on a for large a is somewhat 
intricate, a consequence of the error-function ramp shape. 
For very large values of a — not shown in Fig. [5] — the de- 
cay of |A| 2 becomes slower (the location of the crossover 
depends on the details of the pulse truncation). Although 
an approximate analytic expression for |j4| 2 can be de- 
rived for this large-er limit, the formula is not useful for 
the regimes of interest here. 
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We now turn to the application of the switching error 
formula ((76)) to the CZ gate of Table |U As we have ex- 
plained above in Sec. IIIIE) the initial condition in the 
QVN 4 processor with the minimum fidelity after a CZ 
gate (between qubit q\ and the bus) is the eigenstate 
([ST) . For a perfect CZ gate this state would map to 



-| 100000001), (78) 

but in practice we will find population in other channels 
as well, the dominant error (for an optimal pulse) oc- 
curing in 1 000000002), which has two excitations in the 
bus. This dominant error mechanism involves only a sin- 
gle qubit and resonator, and can be understood entirely 
within the truncated model ©■ In the language of that 
model, where the bare states (interacting eigenstates) are 
written as \qb) (| <?&)), the dominant fidelity loss of the |11) 
eigenstate is caused by nonadiabatic leakage to 1 02) . We 
call this leakage error E2. 

It is possible to understand this dominant |11) —> 1 02) 
leakage error in a further simplified two-channel model 
that only includes the channels |11) and 1 02) . Given the 
strong interaction of the bare |11) state with 1 20) , it is 
not at all obvious that such a simplification is possible. 
However, during most of the switching, the detuning be- 
tween 1 11) and 1 20) is much larger than their interaction 
strength V2gb, so they are effectively decoupled. (And 
while the qubit sits at the w on frequency, the Hamiltonian 
is time-independent.) Therefore a two-channel descrip- 
tion should be possible, altho ugh it will slightly overes- 
timate the excitation of 1 02) . Numerical investigation 
confirms that the 1 20) channel can indeed be disregarded 
except for the level repulsion it produces on the 1 11) state 
(see below) , which is crucial for obtaining an accurate fi- 
delity estimate. 

In the bare {|11), |02)} basis, ([9]) reduces to 

R= f e \Z2SbV ongt (?9) 
Vv2.gb w b / 

At the beginning of the CZ pulse of Table U the qubit 
frequency is 1.0 GHz above the bus. e then decreases to 
uj on w + i] and returns to the detuned configuration 
in the manner of Fig. [3] The total leakage to 1 02) can 
be estimated as twice — because there are two switching 
events, which we assume to contribute incoherently — the 
value of p sw ■ Therefore the error E 2 introduced in (|39|) 
is given by 

E 2 = 2 Psw . (80) 

To evaluate the switching error in this case, we use (|76j) 
with the parameter values 

G = V2g h , (81) 
A on = 77 - V2g h , (82) 
A off = w oH - uj h . (83) 



The value of A on accounts for the level repuls ion c aused 
by the neglected 1 20) state, which causes the |11) chan- 
nel to shift downwards towards 1 02) ; this large effect is 
evident in Fig. [3] [We ignore here a smaller repulsion by 
1 02), which would lead to the addition of a small posi- 
tive correction to (|82p .j A Q ff/27r is always 1.0 GHz for 
the gates listed in Table HI The required \A\ 2 values are 
obtained from Fig. [BJ and are provided in Table U along 
with p sw . 

The minimum fidelity estimate 

F | ( 1 ° 1 s ) t) = l-E 2 = l-2p sw (84) 

is also given in Table U and compared to simulated Q VN 4 
values of F ave and -Fin). We find that (|84|) is a reliable 
predictor of the worst-case fidelity i*ju\ in QVN 4 , con- 
firming that the nonadiabatic switching error is the dom- 
inant fidelity loss mechanism here. Although this error 
will always be present, it can be exponentially suppressed 
by increasing the switching time. 

Finally, we briefly comment on the nonadiabatic 
switching errors between the 1 01) and 1 10) eigenstates, 
which we have argued to be subdominant to the exci- 
tation of 1 02) , but which naively are of the same order. 
There are two reasons why the 1 01) and 1 10) switching 
errors are considerably smaller: First, the matrix ele- 
ment coupling 1 01) and 1 1 0) is a factor of smaller 
than that between |11) and 1 02) , and this factor gets 
squared in (|76[) . And the second — but quantativitely 
more important — reason is that while level repulsion con- 
siderably enhances the 1 02) excitation [recall (|82|) ]. it 
(slightly) suppresses transitions between 1 1) and 1 10) . 

We can estimate the switching errors between the 1 01) 
and 1 10) channels during a CZ gate by using ([76)) with 
parameters 

G = 5 b, (85) 

A on = 7?+^, (86) 
V 

and with A ff/27r = 1 GHz as before. The expression 
([SB") for A on accounts for the level repulsion between 
1 01) and 1 10), which suppresses the switching error, in 
contrast with the strong enhancement indicated in (|82p . 
Considering, for example, the 99.9% CZ gate of Table 
U designed for the rj/2n = 300 MHz qubit, we find that 
\A\ 2 = 2.2 xlO" 3 and 

p sw = 4.4xHT 5 . (87) 

The total error in this case is Ei = 2p sw . T he e stimate 
([57)1 for the switching error between 1 01) and 1 10) is more 
than an order-of-magnitude smaller than that between 
JTTy and ]02), which is p sw = 1.2 x 10~ 3 (see TableO}. 

I. Pulse shape errors 

In Sec. IIII HI we discussed the intrinsic error of the 
qubit-bus CZ gate — assuming an optimal pulse shape — 
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and identified its dominant source as a nonadiabatic 
switching error E 2 . In this section we discuss and quan- 
tify the fidelity loss caused by pulse shape and auxil- 
iary z rotation errors. By a pulse shape error we mean 
that the correct functional form ([lj is seen by the qubit, 
but with values of the parameters t OTL and uj on that de- 
viate from the optimal values. It is possible to develop 
simple analytic models for these error mechanisms (sup- 
ported by numerical simulation) by noting that when the 
fidelity is very close to unity, the different error mecha- 
nisms present contribute independently and can be cal- 
culated separately. We then use these results to estimate 
the experimental pulse-control precision required during 
the implementation of a given 99.9% or 99.99% CZ gate 
to keep any accompanying pulse shape error less than the 
base 10 -3 or 10~ 4 gate error. 

The simplest situation to consider is that where the 
correct values of t on and w on are used, but where the local 
z rotation angles jk {k = 1,2) applied experimentally 
deviate from their optimal values by amounts ipi- -C 1. 
We estimate the resulting fidelity loss by imagining that 
we have achieved a perfect CZ-class gate 



U 



(I 

e 102 



\0 














D j(ei+e 2 ) 



for some phase angles 0^, but then apply z rotation an- 
gles 



Ik = 6fe + 



(89) 



a leading order error E = 1 — F ave given by 



that have errors ipk ■ From ([50)) we find that this leads to 



(90) 



E 



Next we consider i on and w on errors. An error in either 
< on or w on has two consequences, the first is to modify 
the accumulated phases 9^ in (|88|) . and the second is to 
cause population and phase errors on the |11) channel. 
Therefore we consider two types of pulse shape errors, the 
first where i on or w on is changed with no compensating 
changes in the auxiliary z rotation angles, and the second 
where the jk are reoptimized. 

In the first case, in which the error is clearly the largest, 
the resulting error is dominated by the z rotation angle 
error itself, which can be estimated from (|90|) . Chang- 
ing t on by an amount St oa , or w on by an amount 8u) oa , 
changes the accumulated phase of the qubit (recall dis- 
cussion of the qubit reference frame in Sees. 1111 Al and 

Unci) by 

tfl = (Won - ^off) St oa + Su) on ion, (91) 

and that of the bus by 



<p2 = 0. 



(92) 



We note that St on and 8uj on can be positive or negative 
here, and that the total gate time t gate is also (slightly) 
changed by 5t OD . The additional accumulated phase (|9"Tj) 
can be regarded as a rotation angle error because, by 
assumption, it is not compensated by the applied z ro- 
tations (hence the notation). The error angle tf2 is zero 
because of our choice of the local clock/reference frame 
for a resonator. Therefore, an error in either t on or w on 
with no compensating adjustment of the auxiliary z ro- 
tation angles leads to a leading-order fidelity loss of 



E 



(Woff 



■ St. 



2 

on 



t 2 



(93) 



For an order-of-magnitude estimate it is sufficient to ap- 



proximate w on here by w D and t on by t[ 
leading to the simpler estimate 



sudden 



/-/-sudden ^2 



Sui 



2 

on ' 



(94) 



Considering the 99.9% CZ gate of Table Q] designed for the 
i]/2ir = 300 MHz transmon, we estimate from (|94|) that 
a 10 ps error in i on (or a 1 MHz error in w on /27r) would 
lead to a fidelity loss of 7.9x 10~ 4 (or 4.9x 10 -4 ), whereas 
numerical simulation of the same error in QVN4, which 
includes all subdominant processes, yields 4.5 x 10~ 4 (or 
1.2xl0- 3 ). 

Next we consider the case where there is an error in t on 
or w on , but the auxiliary z rotation angles are optimal. 
Here the analysis closely follows that of Sec. IIII B( which 
is based on the truncated qubit-resonator model In 
this situation the fidelity loss is dominated by deviations 
from the ideal evolution 



CZ|11) 



111) 



(95) 



of the 1 11) channel. Pulse shape errors will lead to both 
population and phase errors on the right-hand-side of 
(lS?5"l) . We therefore parameterize the nonideal CZ gate by 



U = phase factor x 



/l 
1 
1 

\0 







— e lS cos 



(96) 



In (|96[) we have assumed perfect auxiliary z rotations 
and have neglected the subdominant errors in the |00), 
1 01) , and 1 10) channels. The population error has been 
written in terms of a rotation angle error 9, which will be 
interpreted (see below) as the deviation from 2ir of the 
rotation angle in the two-dimensional subspace spanned 
by 1 11) and |20). The expression (|M|) does not include 
switching errors because we are evaluating the effect of 
pulse-shape errors on an otherwise perfect, large-tr CZ 
gate. The average fidelity loss E = 1 — F avc associated 
with (|9f)|) is, to leading order, 



E 



20 



16 



(97) 
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What remains is to express the controlled phase error 5 
and rotation error angle 9 in terms of 5t on and <5w on . This 
involves only the two channels |11) and 1 20) of the qubit- 
resonator model ([5]), and we will use the same small-cr 
approximation used in Sec. IIII Bl for our analysis of the 
|11) channel dynamics. In the {|11), 1 20) } basis, © can 
be written as 



\V2g h 2e-r) 



The eigenstates of (|98| are 



|11) = cos| |11) -sin| |20), 
py = cos^|20)+sin^|ll), 



where 



£ = arctan 

and the energies are 

E n = e + oj h + ■ 



V £ - Won / ' 



(98) 

(99) 
(100) 

(101) 



e - uj 



(V2g h ) 



E 2 o = e + uj h 



(102) 



(103) 



Here we have used the expression (fTTj) for ui on , which is 
appropriate for the model (|9"5|) . 

The analysis below follows that of the |11) channel 
evolution given Sec. IIII Bt except here we introduce a 
timing error 6t on and a tuning error <5u; on . Starting in 




the strongly detuned configuration with e 



> gfc 



in the eigenstate |11) ~ |11), and quickly switching to 
e = uj on + <5w on , leaves the system in the state 



|11) = cos ^ |11) 



sin%|20). 



(104) 



Here £ n = arctan(2\/2gb/<5k-'on), and the eigenstates in 
(|104| are for e = uj on + 5uj on . Note that (|104[) reduces 
limit. Evolution with e fixed at 



to 



in the Su> n 



2n 



20 



V / (<5w on /2) 2 + (y2, 



(105) 



2.9b 



would implement a 2ir rotation in the {|11),|20)} sub- 
space, returning to |11) with a phase shift that depends 
on Suj on . We intentionally introduce a t on pulse shape 
error and instead evolve for a time 



t = 



St a 



(106) 



after which we detune and find the final state 

- e lS cos - 1 11) + phase factor x \/Eq 1 20) 



e iS cos- 1 11) + phase factor x ifEi,\20), (107) 



where 



and 



6 = - 



■Suj„ 



2V2g h 



En 



(108) 



(109) 



These expressions are valid to leading order in 5w on or 
6t on , neglecti ng cr oss terms. Here E# is the probability 



of leakage to 1 20) resulting from a t on error, which, as 
discussed above, causes a rotation error of angle 9. An 
alternative estimate for 9 (and hence Kg) is 



St, 



x 2tt, 



(110) 



which also gives (|109[) [after using the sudden limit result 
(HI for t on ]. Note that the leakage error E# is indepen- 
dent of 6uj on (to this order), enabling the phase 8 to be 
intentionally adjusted by varying w on only. Doing this 
generates (approximately) gates of the form 
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\0 e^+V 
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(111) 



for small S, with alternative controlled phases [46j . Gates 
(approximately) of the form (|1 1 1|) with arbitrary — not 
necessarily small — values of S can be implemented by 
varying both u> on and t on away from their optimal val- 
ues. 

Referring again to the 99.9% CZ gate of Table Q] de- 
signed for the r]/2ir = 300 MHz transmon, we estimate 
from ([571), (fTUSj) . and (fTM)) that a 10 ps error in t on 
(or a 1 MHz error in w on /27r) with optimal z rotations 
would lead to a fidelity loss of about 4.0 x 10~ 6 (or 
9.1x 10 -5 ), whereas numerical simulation of the same er- 
ror in QVN4, which includes all subdominant processes, 
yields 3. Ox 10~ 6 (or 1.1 xl0~ 4 ). 

Finally, it is interesting to use the above estimates to 
bound the magnitude of the allowable t on and w on errors, 
such that the resulting pulse shape errors ([97)) are sub- 
dominant to the 10 -3 or 10~ 4 base gate error. This is 
done in Table HTT1 For example, a 99.9% CZ gate from Ta- 
ble[T]with a t on error of 160 ps will have an additional in- 
trinsic error of 10 -3 (and a total error of about 2 x 10~ 3 ). 
Current experimental limitations on the control of t on 
and uj on are considerably better than that required to 
suppress pulse shape errors below the 10 -4 level. 
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TABLE III. CZ pulse shape precision requirements. The 
bounds listed in the t on column assume that this is the only 
type of pulse parameter inaccuracy present, with an (esti- 
mated) error given in the first column, and that the auxiliary 
z rotation angles are reoptimized and implemented perfectly. 
The Won bounds are defined analogously. The error E is de- 
fined in {§3). 



E 


t on precision 


Won precision 


icr 3 

10" 4 


160 ps 
50 ps 


4 MHz 
1MHz 



J. CZ gate with 99.999% fidelity 

Higher fidelities are also possible with the pulse shape 
([T]). An example is provided in Table ED for the 300 MHz 
qubit. For this design we did not perform a separate <?b 
optimization for this higher fidelity, but instead used the 
value from Table U optimized for the lower fidelities. 



TABLE IV. Optimal QVN 4 gate fidelity for a Strauch CZ gate 
between qubit q\ and the bus. 
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FIG. 6. (Color online) QVN 4 mode diagram for the CZ23 
gate. Gaussian filtering of the pulse is not shown. 



IV. ADDITIONAL CZ GATES 

A. CZ between bus and qubit 

The main focus of this paper is the CZ gate between 
qubit q\ and the bus in the QVN 4 processor. Results 
for the other qubits are very similar, with the worst case 
being 54, because the detuning to memory during the 
gate is slightly smaller. We find that the intrinsic fidelity 
of the 99.9% CZ gate for the 300 MHz qubit given in 
Table U changes from 99.928% to 99.925% if qubit g 4 is 
used instead of q±. 



B. CZ between two qubits in QVN 4 

In this section we explain how to perform a CZ gate 
between two qubits — or more precisely, between two 
memories — in the QVN architecture. Such an operation 
is not elementary, as it can be composed of the qubit-bus 
CZ combined with MOVE gates. (There are also pro- 
posals for the direct implementation of a qubit-qubit CZ 
gate in a QVN device Figure [6] shows the ex- 

perimental protocol for implementing the gate CZ23 be- 
tween qubits (72 and (73, suppressing auxiliary z rotations, 
and with all data starting and ending in memory. This 
mode diagram shows the time-dependence of all 9 device 
frequencies. The color indicates whether the qubit or res- 
onator would be in the ground state (blue) or possibly 
excited state (red) in the weakly coupled limit. Modes 



are colored red if there is a finite occupation of the |1) 
state (in the weakly coupled limit), for some choice of 
initial conditions. All qubits are initially parked at the 
strongly detuned frequency of 10 GHz. Horizontal lines 
1-4 represent memories, and b is the bus. No red or red 
and blue lines with first-order or second-order couplings 
cross (to avoid Landau-Zener transitions), and no more 
than one qubit is occupied at any time (to avoid second- 
order qubit-qubit interactions mediated by the bus). 

Beginning with the four memory registers in an arbi- 
trary (possibly entangled) state, the bus is loaded by a 
5 ns MOVE gate from 7713 —> followed by a 10 ns MOVE 
to the bus. These are approximate gate times (time esti- 
mates for these gates are given in Sec. |V] and a concrete 
example is provided below) . Qubit q 2 is then loaded and 
tuned to the frequency w on determined by optimization. 
This central portion of the gate is close (but not exactly 
the same as) the qubit-bus CZ gate of Table U which, for 
the purposes of comparison, is shown in Fig. [7] 

We simulated the gate CZ23 shown in Fig. |H1 using the 
99.99% CZ of Table [D for a 300 MHz qubit, and starting 
with the memory register in the GHZ state 



|0000) + |1111) 
V2 ' 



(112) 



The MOVE gates also have fidelities around 99.99%. Note 
that -ffjdie and the associated computational basis states 
(interacting eigenfunctions of -ffidic) ar< 3 different at the 
beginning and end of each MOVE gate. The total CZ23 
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FIG. 7. (Color online) CZ gate between qi and bus. This 
diagram describes the gate of Table |I] 



gate time is 



t 



gate 



55 ns, 



(113) 



and the final state fidelity (overlap squared) is 

F = 99.94%. (114) 

In addition to the 23 ns qubit-bus CZ gate, there are four 
MOVE to/from memory operations, each taking about 
3.5 ns, and two MOVE to/from bus gates, each taking 
about 9 ns. There are also local z rotations (not shown) 
between each gate. 



A few remarks about the encouraging result (jl 14[) arc 
in order: The seven elementary gates making up the CZ23 
operation are optimized individually to an error of about 
10~ 4 , and then combined without any additional opti- 
mization of the composite pulse sequence or control pa- 
rameters, respecting the modularity required by scalable, 
gate-based quantum computation. And the total intrin- 
sic error E = 1 — F = 6 x 1CP 4 implied by (|114j) is con- 
sistent with a linear (incoherent) accumulation of errors 
with number of elementary steps Ej nc = 7 x 10 -4 , but is 
not consistent with a quadratic (coherent) accumulation 
E coh = 7 2 x 1(T 4 . 



C. Beyond QVN: CZ gate between directly coupled 
superconducting qubits 

The main focus of this work is the CZ gate between 
a directly coupled qubit and resonator in a multi-qubit 
QVN device. However many of our results — especially 



the error analysis — will also be relevant for hardware 
designs incorporating pairs of directly coupled super- 
conducting qubits, the system originally considered by 
Strauch et al. 42j . Here we summarize the principal dif- 
ferences between the qubit-bus CZ gate of Sec. IIIII and 
the directly coupled qubit-qubit gate. For the latter case 
we assume qubits with anharmonic detunings 771 and 772 
[denned as in ([7])] and a purely transverse (off-diagonal) 
capacitive coupling with interaction strength g. 

The first difference concerns a frequency asymmetry of 
the qubit-bus gate. Recall from Sec. IIIIBI that the CZ 
gate is implemented by decreasing the qubit frequency e 
from a value far above Wb to (approximately) 



e = ui h + ?y, 



(115) 



where 77 > is the qubit anharmonicity. By contrast, 
the qubit-qubit CZ gate can be implemented either by 
decreasing the frequency ei of qubit 1 from a value far 
above that of a second qubit with a (fixed) frequency €2, 
until 



ei = £2 



(116) 



which is directly analogous to (|115p . or from below by 
increasing to 



ei = £2 



(117) 



The conditions (|116[) and (|117[) specify the crossings of 
the bare |11) state with 1 20) and |02), respectively (in 
the basis |<7i<72))- The frequency asymmetry of the qubit- 
bus gate is a consequence of the harmonic spectrum of 
the bus and can be understood from Fig. [3J which shows 
that when |11) reaches the 1 02) crossing from below, 1 01) 
and 1 10) also become degenerate [as expected from (|117[) 
when 772 — > 0] . This would result in unwanted phase shifts 
of the 1 01) and 1 10) channels, as well as large switching 
errors between them. 

The second major difference between the qubit-bus 
and qubit-qubit gates is that the additional anharmonic- 
ity in the latter case further suppresses the nonadia- 
batic switching errors and leads to better gate perfor- 
mance. This can be understood from the analysis of 
Sec. lIIIHl noting that in the qubit-qubit case, adiabatic- 
ity of the 1 1 1 ) channel is protected by an energy gap of 
size 771 + 7/2 — V%g, where g is the qubit-qubit interaction 
strength. Let's estimate the 1 1 1 ) — > 1 02) switching er- 
ror E2 for a qubit-qubit CZ gate at the upper frequency 
([TT7J1) . with 771 = 772 = 2tt x 300 MHz, g = 45 MHz, and 
t, 

parameter values 



ramp = 7 ns. In this application we use formula ([76]) with 



G = V2g, 



771 + 772 - V2g, 
2tt x 1 GHz. 



(118) 
(119) 
(120) 



We note from (| 1 19[) that the anharmonicity suppressing 
the 1 1 1) — > 1 02) switching error is effectively doubled in 
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the qubit-qubit system. With these parameters we obtain 
\A\ 2 = 5.8 x 1CT 6 and p sw = 8.2 x 1CT 8 . 

Although the |11) — > 1 02) switching error is greatly 
reduced in the qubit-qubit CZ gate, the actual gate fi- 
delity does not fully benefit from this reduction. This is 
because the dominant intrinsic error in the qubit-qubit 
gate is the switching error Ei between |01) and |10), or 
the reverse, which is subdominant in the qubit-bus case 
(see Sec. IrnHl) . In fact, the |01) O jlO) switching error 
estimate (|87|) also applies to the qubit-qubit system (with 
r h =t] 2 = 2?rx300MHz, g = 45 MHz, and i ramp = 7ns), 
resulting in an estimated minumim fidelity of 



F^ ) = l-E 1 = l-2 Pa 



= 99.991%. 



(121) 



the (approximate two-parameter) MOVE to/from mem- 
ory gate becomes significant (this is determined by the 
qubit-bus detuning because the dominant error is nona- 
diabatic leakage to the bus), whereas if uj g is too high 
the fidelity of the qubit-bus CZ degrades (because de/dt 
increases). We find that 7.5 GHz works well. At least 
500 MHz is required between w ff/27r and the memory 
band to keep the MOVE errors (to/from the bus) under 
control. Thus we arrive at the memory frequencies given 
inTable|Hl 

Having obtained prospective resonator frequencies, we 
turn to couplings. The most frequently used gate is ex- 
pected to be the MOVE to/from memory, which must be 
as fast as possible. Figure [B] includes four examples. The 
gate time for this operation is approximately 



V. SYSTEM OPTIMIZATION 

In this section we discuss an approach for choosing 
optimal QVN n device parameters. This is a complex 
global optimization problem that we will solve in a simple 
but approximate way, emphasizing the main ideas of the 
procedure instead of its most precise implementation. 

First we consider resonator frequencies. The QVN„ 
processor includes n memory resonators, with frequencies 
w m i,w m 2, • • • . These need to be mutually detuned (to 
lift degeneracies), and for simplicity we space them by 
100 MHz (a smaller value could be used). The band of 
memory frequencies itself needs to be well detuned from 
the bus to keep the idle error (to be discussed below) in 
check. 

Because the qubit frequency during a qubit-bus CZ 
gate must approach the bus frequency from above 
(Sec. HVCj) . the bus frequency must be below the mem- 
ory band. The choice of bus frequency therefore de- 
termines the lowest transition frequency that needs to 
be accessible by a qubit. Specifically, the qubits will 
need to tune 500 MHz or so below the bus (see Fig. [B]) . 
However, the minimum transition frequency may be con- 
strained by qubit design (in addition to other consider- 
ations). In the tunable- Ej transmon, for example, this 
minimum frequency depends on the qubit anharmonicity 
77. Here we will choose a minimum qubit frequency and 
corresponding bus frequency appropriate for a 300 MHz 
transmon. This leads to our choice of 6.5 GHz for the 
bus frequency. Optimal resonator frequencies for smaller 
77 are unchanged, whereas for larger 77 they need to be 
rigidly shifted upward in frequency. In particular, sys- 
tem frequencies for a 400 MHz transmon will be shifted 
upward in frequency by about 2 GHz. Apart from this 
large but simple change, we expect the system optimiza- 
tion results, such as <?b values, to be valid for the 400 MHz 
case as well. 

The frequency to g can be viewed as defining a bound- 
ary between MOVE and CZ gates, or between consecutive 
MOVE gates (see Fig. [BJ for example). It is also natural 
to perform single-qubit operations with microwave pulses 
at the qubit frequency u> g . If lo r is too low, the error of 



+ <ramp + 1 nS. 



(122) 



The first term is the 7r rotation time, and the sec- 
ond and third are switching times (the detuning ramp 
can be fast because the qubit is unoccupied). Choos- 
ing g m /2ir = 100 MHz makes the first term 2.5 ns. It 
might be possible to increase g m further, but suppresss- 
ing the resulting idling error (see below) would require an 
even higher empty-qubit parking frequency. The value 
of iramp is determined by the desired MOVE gate fi- 
delity. Because the dominant error is nonadiabatic ex- 
citation of the bus, we can estimate it using our ex- 
pression (|76|) for the switching error p SVJ , with G = <?b, 
A on /2n = 1.0 GHz, and A off /27r = 1.8 GHz (correspond- 
ing to 7714, the worst case). These values depend on our 
initial resonator frequency assignments. From Fig. [5] we 
obtain \A\ 2 = 0.03 (9.8 x 10" 4 ) for a Ins (2 ns) ramp. 
Considering the largest (worst case) value for <?b/27r of 
60 MHz gives p sw = l.lxl0~ 4 (3.5xl0~ 6 ) for a I ns (2 ns) 
ramp. Thus we conclude that the MOVE to/from mem- 
ory can be done in about 5 ns if g m /2n = 100 MHz. (Here 
we assumed the simplest 2+1-parameter MOVE to/from 
memory gate, having two pulse-shape parameters and 
one auxiliary z rotation angle. It is also possible to im- 
plement this gate with even higher fidelity with 4+1 pa- 
rameters [34j •) 

The bus coupling is found by the following "g- 
optimization" procedure: Consider the set of discretized 
<7b/27r values, varying from 10 to 100 MHz in steps of 
1 MHz. For each value of g\>, calculate the minimum value 
of tramp and the associated £ ga t c required to achieve a tar- 
get fidelity, say 99.9%. We do this by stepping through 



L ramp 



values, estimating the fidelity using (1761) and 



from Sec. IIIIH1 which is very efficient, then confirming 
through a full optimization on QVN-p We then obtain, 
for each g\>, the gate time of a 99.9% CZ gate, or equiva- 
lently, the function 



t 



(99.9%) 

gate 



(123) 



The function (j 123[) gives the time required for a CZ gate 
with a given target fidelity as a function of Strict g- 
optimization requires choosing g^ to minimize i ga te, and 
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this procedure leads to the best performance in any given 
situation. However, carrying this out leads to a different 
gh for each target fidelity. Fortunately, the curvature at 
the minimum in (|123p is small and the simpler <?b values 
reported in Table HI which are independent of the target 
fidelity and are also constrained to be multiples of 5 MHz, 
lead to very little performance loss, typically 1-2 ns in 
gate time. 

Having obtained prospective resonator frequencies and 
couplings, we choose the empty qubit parking frequency 
Wpark to control the idle error 

E = (f! zz i) V, (124) 

where £lzz is the effective a z ® a z coupling frequency 
between a memory resonator and the bus, induced by 
their shared qubit [34]. The n-dependence in (|124p as- 
sumes the worst case. Assuming w pa rk/27r = 10 GHz, 
g b /27T = 60 MHz (the worst case), and rj/2ir = 400 MHz 
(also the worst case) leads to 

— = -0.881 kHz. (125) 

It will be necessary to keep (|124j) less than the fault- 
tolerant threshold during a (potentially) long error cor- 
rection cycle. If we assume t — 1 /js, the idle error in 
QVN 4 is 4.9 x 10~ 4 , which is acceptable. Reducing the 
parked qubit frequency to 9.5 GHz (9.0 GHz) increases 
the idle error to 3.8 x 10~ 3 (7.3 x 10~ 2 ). 

We are now able to calculate the gate time of the 
MOVE to/from bus operation. The gate time is approx- 
imately 

7T 

■Z h £ramp + 1 nS, (126) 

25b 

where the first term is between 4 and 8 ns for the bus 
couplings in Table HI As before, t ramp is determined by 
the desired gate fidelity. The dominant error is nona- 
diabatic transition to memory, which we estimate us- 
ing (|76|) with G = g m , A on /2n — 0.5 GHz (the worst 
case), and A s/2tt = 1.5 GHz. Note that this error is en- 
hanced by the large value of g m . From Fig. [5] we obtain 
\A\ 2 = 5.9xl0~ 2 (1.6xl0~ 2 ) for a 2 ns (3 ns) ramp. Then 
p sw =2 Ax 10~ 4 (6.4 x 10" 4 ) for a 2 ns (3 ns) ramp. Thus 
we conclude that the MOVE to/from bus takes between 
7 and 12 ns, depending on the actual value of g^ and on 
the desired fidelity. 

Finally, we confirm that the assumed qubit parame- 
ters are compatible with transmons. In the large Ej/Eq 
transmon regime, the qubit frequency e and anharmonic- 
ity 77 are given by [|| 

e = ^8EjE c - E c (127) 

and 

r/ = E c . (128) 



We assume a split-junction flux-biased Cooper-pair box 
so that Ej is tunable [43j . (Note that the tunablc- 
Ej transmon is sensitive to flux noise, which will de- 
grade T 2 .) Combining (|127p and (I128[) leads to the re- 
lation e = r)(y/8Ej/E c - 1) plotted in Fig. ® for 300 
and 400 MHz anharmonicity. Because Ej / Eq needs to 
be above about 50 to effectively suppress charge noise, 
we see that the 300 MHz transmon can have a transi- 
tion frequency as small as 5.5 GHz, whereas the 400 MHz 
transmon has a minimum frequency of about 7.5 GHz. 
Our choice of bus frequency is indeed consistent with the 
300 MHz transmon, whereas Wb/27r and the entire spec- 
trum of device frequencies would have to be increased by 
about 2 GHz for the 400 MHz transmon. 




FIG. 8. (Color online) Transition frequency for transmon with 
300 and 400 MHz anharmonicity. 



VI. CONCLUSIONS 

We have investigated the problem of CZ gate design for 
the superconducting QVN architecture based on a real- 
istic, two-parameter filtered rectangular pulse (fTJ). The 
gate operates by using the non-computational qubit |2) 
state as first proposed by Strauch et al. [42|]. The main 
results are summarized in Table HI The use of interact- 
ing eigenfunctions as computational basis states, and the 
use of auxiliary z rotations on the qubits and bus are 
critical to obtaining this high performance. The ability 
to perform SU(2) operations on the bus resonator is non- 
trivial because they cannot be implemented directly with 
microwave pulses or frequency excursions; instead, they 
must be compiled into future qubit rotations. 

Our investigation is limited to and constrained by the 
pulse profile (TTJ). Fidelity optimization determines the 
amount of time to spend switching the qubit and how 
long to stay at the "on" frequency. We find that this 
pulse shape correctly captures the relevant pulse degrees- 
of-freedom for fidelities up to about 99.99%. For higher 
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fidelity gates, which are slower and which would require 
larger coupling, alternative few-parameter pulse shapes 
where most of the time is spent switching, could become 
competitive with (Q~]). Although we are able to obtain 
fidelities well above 99.99% with profile ((T|), it is pos- 
sible that alternative few-parameter pulse shapes would 
be able to achieve the same intrinsic fidelity in less gate 
time. 

One can also consider more complex pulse shapes with 
many control parameters, which can achieve nearly per- 
fect intrinsic fidelity in a time i ga t (depending on gh) 
significantly shorter than obtained with pulse shape ((T|) . 
Egger et al. [5!| have recently investigated this optimal 
control approach (using the gradient pulse shape engi- 
neering method of Khaneja et al. [60]), and have obtained 
about a factor of two speedup for a qubit-qubit CZ gate 
similar to that of Sec. IIVB1 This approach clearly war- 



rants further investigation and experimental implemen- 
tation. 

We did not include the effects of decoherence (or flux 
noise) in this work. However, an order-of-magnitude es- 
timate of the Ti decay error E « t gate /Ti suggests that it 
should be possible to demonstrate a 99.9% CZ gate with 
existing transmon qubits, which would be an important 
step towards the development of fault-tolerant quantum 
computation. 
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